Differential haplotype expression in class I MHC genes during SARS-CoV-2 infection of human lung cell lines.

Introduction: Cell entry of SARS-CoV-2 causes genome-wide disruption of the transcriptional profiles of genes and biological pathways involved in the pathogenesis of COVID-19. Expression allelic imbalance is characterized by a deviation from the Mendelian expected 1:1 expression ratio and is an important source of allele-specific heterogeneity. Expression allelic imbalance can be measured by allele-specific expression analysis (ASE) across heterozygous informative expressed single nucleotide variants (eSNVs). ASE reflects many regulatory biological phenomena that can be assessed by combining genome and transcriptome information. ASE contributes to the interindividual variability associated with the disease. We aim to estimate the transcriptome-wide impact of SARS-CoV-2 infection by analyzing eSNVs. Methods: We compared ASE profiles in the human lung cell lines Calu-3, A459, and H522 before and after infection with SARS-CoV-2 using RNA-Seq experiments. Results: We identified 34 differential ASE (DASE) sites in 13 genes (HLA-A, HLA-B, HLA-C, BRD2, EHD2, GFM2, GSPT1, HAVCR1, MAT2A, NQO2, SUPT6H, TNFRSF11A, UMPS), all of which are enriched in protein binding functions and play a role in COVID-19. Most DASE sites were assigned to the MHC class I locus and were predominantly upregulated upon infection. DASE sites in the MHC class I locus also occur in iPSC-derived airway epithelium basal cells infected with SARS-CoV-2. Using an RNA-Seq haplotype reconstruction approach, we found DASE sites and adjacent eSNVs in phase (i.e., predicted on the same DNA strand), demonstrating differential haplotype expression upon infection. We found a bias towards the expression of the HLA alleles with a higher binding affinity to SARS-CoV-2 epitopes. Discussion: Independent of gene expression compensation, SARS-CoV-2 infection of human lung cell lines induces transcriptional allelic switching at the MHC loci. This suggests a response mechanism to SARS-CoV-2 infection that swaps HLA alleles with poor epitope binding affinity, an expectation supported by publicly available proteome data.

Genomic Perspectives on the Emerging SARS-CoV-2 Omicron Variant.

A new variant of concern for SARS-CoV-2, Omicron (B.1.1.529), was designated by the World Health Organization on November 26, 2021. This study analyzed the viral genome sequencing data of 108 samples collected from patients infected with Omicron. First, we found that the enrichment efficiency of viral nucleic acids was reduced due to mutations in the region where the primers anneal to. Second, the Omicron variant possesses an excessive number of mutations compared to other variants circulating at the same time (median: 62 vs. 45), especially in the Spike gene. Mutations in the Spike gene confer alterations in 32 amino acid residues, more than those observed in other SARS-CoV-2 variants. Moreover, a large number of nonsynonymous mutations occur in the codons for the amino acid residues located on the surface of the Spike protein, which could potentially affect the replication, infectivity, and antigenicity of SARS-CoV-2. Third, there are 53 mutations between the Omicron variant and its closest sequences available in public databases. Many of these mutations were rarely observed in public databases and had a low mutation rate. In addition, the linkage disequilibrium between these mutations was low, with a limited number of mutations concurrently observed in the same genome, suggesting that the Omicron variant would be in a different evolutionary branch from the currently prevalent variants. To improve our ability to detect and track the source of new variants rapidly, it is imperative to further strengthen genomic surveillance and data sharing globally in a timely manner.

Host genetic susceptibility underlying SARS-CoV-2-associated Multisystem Inflammatory Syndrome in Brazilian Children.

BACKGROUND: Multisystem Inflammatory Syndrome in Children (MIS-C) is a life-threatening complication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, which manifests as a hyper inflammatory process with multiorgan involvement in predominantly healthy children in the weeks following mild or asymptomatic coronavirus disease 2019 (COVID-19). However, host monogenic predisposing factors to MIS-C remain elusive. METHODS: Herein, we used whole exome sequencing (WES) on 16 MIS-C Brazilian patients to identify single nucleotide/InDels variants as predisposition factors associated with MIS-C. RESULTS: We identified ten very rare variants in eight genes (FREM1, MPO, POLG, C6, C9, ABCA4, ABCC6, and BSCL2) as the most promising candidates to be related to a higher risk of MIS-C development. These variants may propitiate a less effective immune response to infection or trigger the inflammatory response or yet a delayed hyperimmune response to SARS-CoV-2. Protein-Protein Interactions (PPIs) among the products of the mutated genes revealed an integrated network, enriched for immune and inflammatory response mechanisms with some of the direct partners representing gene products previously associated with MIS-C and Kawasaki disease (KD). In addition, the PPIs direct partners are also enriched for COVID-19-related gene sets. HLA alleles prediction from WES data allowed the identification of at least one risk allele in 100% of the MIS-C patients. CONCLUSIONS: This study is the first to explore host MIS-C-associated variants in a Latin American admixed population. Besides expanding the spectrum of MIS-C-associated variants, our findings highlight the relevance of using WES for characterising the genetic interindividual variability associated with COVID-19 complications and ratify the presence of overlapping/convergent mechanisms among MIS-C, KD and COVID-19, crucial for future therapeutic management.

Identification and characterisation of SARS-CoV-2 and Human alphaherpesvirus 1 from a productive coinfection in a fatal COVID-19 case.

BACKGROUND: During routine Coronavirus disease 2019 (COVID-19) diagnosis, an unusually high viral load was detected by reverse transcription real-time polymerase chain reaction (RT-qPCR) in a nasopharyngeal swab sample collected from a patient with respiratory and neurological symptoms who rapidly succumbed to the disease. Therefore we sought to characterise the infection. OBJECTIVES: We aimed to determine and characterise the etiological agent responsible for the poor outcome. METHODS: Classical virological methods, such as plaque assay and plaque reduction neutralisation test combined with amplicon-based sequencing, as well as a viral metagenomic approach, were performed to characterise the etiological agents of the infection. FINDINGS: Plaque assay revealed two distinct plaque phenotypes, suggesting either the presence of two severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strains or a productive coinfection of two different species of virus. Amplicon-based sequencing did not support the presence of any SARS-CoV-2 genetic variants that would explain the high viral load and suggested the presence of a single SARS-CoV-2 strain. Nonetheless, the viral metagenomic analysis revealed that Coronaviridae and Herpesviridae were the predominant virus families within the sample. This finding was confirmed by a plaque reduction neutralisation test and PCR. MAIN CONCLUSIONS: We characterised a productive coinfection of SARS-CoV-2 and Herpes simplex virus 1 (HSV-1) in a patient with severe symptoms that succumbed to the disease. Although we cannot establish the causal relationship between the coinfection and the severity of the clinical case, this work serves as a warning for future studies focused on the interplay between SARS-CoV-2 and HSV-1 coinfection and COVID-19 severity.